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ABSTRACT 



Context. Vibrationally excited emission from the SiO and H2O molecules probes the innermost circumstellar envelopes of oxygen-rich red 

giant and supergiant stars. VY CMa is the most prolific known emission source in these molecules. 

Aims. Observations were made to search for rotational lines in the lowest vibrationally excited state of HtO. 

Methods. The APEX telescope was used for observations of H2O lines at frequencies around 300 GHz. 

Results. Two vibrationally excited H9O lines were detected, a third one could not be found. In one of the lines we find evidence for weak 
maser action, similar to known (sub)millimeter v 2 =l lines. We find that the other line's intensity is consistent with thermal excitation by the 
circumstellar infrared radiation field. Several SiO lines were detected together with the H 2 lines. 
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1. Introduction 

In the past, millimeter-wavelength rotational lines from 
the lowest vibrationally excited state of water (H2O), 
the V2=l bending mode, have been detected from the 
inner most envelopes of oxyg en-rich giant and supergiant 
stars. lMenten & Melnickllll989l) used the IRAM 30m telescope 
to discover the V2=l, Jk„k<.=550 — 643 line near 232GHz 
toward the high mass-loss supergiant VY CMa and the nearby 
semiregular variable WHya. Toward the former they also 
detected the 96 GHz 440 - 533 transition. Negative results were 
obtained for a number of other stars. All of the observed 
spectra contained a narrow (width as 1 km s _1 ) feature, clearly 
indicating maser action. This feature also shows up in the 
spectrum of the 658 GHz V2=l, 1 1 — lot ortho rotational 
ground-state transition, which has been found to show very 
strong maser action not only in VY CMa and WHya, but in 
a variety of other (super)giant stars jMenten & Young|ll995l) . 
The excitation of these masers is discussed in Sect. 14. II 

There is also the interesting possibility that thermal 
emission from hot water has been detected: The 232 GHz 
spectra of, both, VY CMa and W Hya show, in addition to the 
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maser spike, broad emission over the velocity range covered by 
the circumstellar outflows of these objects. 

The target of the observations is the peculiar red 
supergiant VY Can is Majoris. With a luminosity of 3.5 10 5 L Q 
(Sopka et al. 1985) and a mass- loss rate of 2-3 1O~ 4 M yr" 1 
(Dan chi et alJ Tl994) VYCMa is a remarkable object by any 
standard. For our discussions it is relevant to note that a 
plethora of maser lines from (ground-state and V2=l) water 
and all the Si isotopes of silicon monoxide (from the v = to 
5 states) have been detected from this object. At far-infrared 
wavelengths its spectrum is completely dominated by H2O 
lines including several lines from the vz=l state jNeufeld et alJ 
Il999h . 

We used APEX to search for the H2O V2 = l, 6(,\ —1st, 
660-753, and 523-616 lines near 294, 297, and 336 GHz. 
Our observations are summarized in Sect. |2 In Sect. [3] we 
report the detection of the first and the third, but not of the 
second transition and present the observed spectra along with 
spectra of simultaneously observed silicon monoxide lines. 
After summarizing (in Sect. 14. ll > our phenomenological picture 
of vibrationally excited H2O maser excitation, in Sect. l4.2l we 
investigate the possibility whether either or both of the detected 
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lines may be the result of thermal excitation and come to the 
conclusion that this could be the case for one of them. 

2. APEX observations and data reduction 

Our observations were made in 2005 July/August under 
generally excellent weather conditions with the 12 m Atacama 
Pathfinder Experiment telescope (APEX 1 ; see Giisten et al.; 
this issue). 

The APEX 2a facility receiver was used (Risacher et al., 
this issue), which covers the 278-375 GHz frequency range 
that contains the frequencies of several H2O lines including 
the ones we observed. Two isotopic silicon monoxide lines 
were in the same band as the 336.2 GHz line. A number 
of other lines were also observed including the 29 SiO v=0, 
7=8-7 line. All lines considered here are listed in Tableland 
are discussed in Sect. [3] Calibration was obtained using the 
standard chopper wheel technique. The radiation was analyzed 
with the MPIfR Fast Fourier Transform spectrometer that 
provides 16384 frequency channels (Klein et al., this issue) 
over the 1 GHz intermediate frequency bandwidth. To increase 
the signal to noise ratio, the spectra were smoothed to effective 
velocity resolutions appropriate for the measured linewidths, 
typically 1 km s _1 . To check the telescope pointing the receiver 
was tuned to the 345.8 GHz CO J =3-2 line. We measured 
5 point rasters in azimuth and elevation centered on VY CMa 
itself and pointing corrections were determined by least squares 
fits. The pointing was found to be accurate to within ~ 4", 
acceptable given the beam size, 6*b, which is 18" FWHM at 
345.8 GHz. We present our line intensities in a main-beam 
brightness temperature (Tmb) scale and assumed a main-beam 
efficiency, 77MB, equal to the observationally determined value 
at 345 GHz, i.e. 0.70. We estimate that our intensity scale is 
accurate to within 15%. The FWHM beam size, #b, is 19" at 
336 GHz and 21" at 294 GHz. 

3. Results 

Fig. H presents our observed spectra; see also Table ^ We 
detect the two ortho-H20 lines (523 — 616 and 6gi - 752), but 
not the latter line's para H20-equivalent (660 - 753). Given that 
the 293 and 297 GHz lines Einstein A-values and energies 
above ground are very similar, we would expect an intensity 
ratio of »3:1 for these lines if they were thermally excited. 
For the integrated intensities we obtain a lower limit of 8.2 
for the ratio of the ortho to the para line, suggesting weak 
maser action that boosts the ortho line's intensity more than 
the para line's. The silicon monoxide lines are all between 
the /=8 and 7 rotational levels but within different vibrational 
states of vastly different energy above ground state. The 3l) SiO, 
v=l line clearly shows maser action expressed in multiple 
emission spikes while all the other lines' profiles are consistent 
with (quasi)thermal emission. All lines are centered around 

1 This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 



20 km s which is within a few km s of the stellar velocity 
(whose exact value has an uncertainty of a few km s _1 ). The 
velocity range covered by the 29 SiO, v=0 line is widest and is 
comparable to that observed in ot her thermally excite d lines 
of SiO and other molecules (e.g. iNvman et alJ Jl99 2)). The 
v=3 29 SiO line and the two detected H2O lines are reasonably 
well fit by Gaussians. The very narrow width of v=3 29 SiO 
line (~ 8 km s _1 ) indicates maser action. The H2O lines have 
about half the width of the 29 SiO, v=0 line (and much poorer 
signal-to-noise ratio) and their shapes do not readily suggest 
masing, but, possibly, that they are formed in the innermost 
envelope where the stellar outflow has not yet reached its 
terminal velocity. 

4. Discussion 

4.1. Vibrationally excited water masers 

lAlcolea & Mentenl 

lll993l) noted that the V2=l maser lines described above seem 
to result from a systematic overpopulation of (what they call) 
the "transposed backbone" levels of the energy level diagram, 
(see their Fig. 1), i.e. levels with quantum numbers K a =J and 
K c =0 and 1 for para- and ortho-H20, respectively, with respect 
to their neighboring levels. This situation is just opposite from 
that found in the ground state, where many of the "backbone" 
levels (Ka=0 or 1 and K C =J) may become overpopulated. 
This systematic inversion may be explainable by means of 
differential trapping in the radiative vibrational decays of the 
levels involved in those maser transitions and, possibly, by 
pumping to and decay from higher vibrational states. Detailed 
calculations to support this are still extant and not yet feasible, 
since collisional rates for vibrational transitions are unknown 
and can only be guessed (see Sect. l4.2l >. 

The mechanism described above leads to the excitation of 
the known 440 - 533, 5so - 643, and 1 10 - lot masers at 96, 232, 
and 658 GHz, respectively, and predicts inversion also, a.o., of 
the 294 and 297 GHz 6ei - 7s2 and 6eo - 753 lines, of which we 
only detect the former one. However, given the uncertainties in 
our understanding of the maser excitation process it is by no 
means clear that these lines should be inverted at all. Certainly, 
there is no a priori reason to expect maser emission in the 
336 GHz line. 

Given the spectral similarities between SiO maser lines 
and V2=l H2O lines pointed out bv lMenten & Melnickl Jl99ll) 
it is highly likely that V2=l H2O maser emission arises 
from the same region as the multi-isotope/multi-v/multi-/ SiO 
maser emission modeled by iGonzalez-Alfonso & Cernicharol 

GUI 2 . 

It may thus be expected that for some of the known 
V2=l lines part or all the emission is masing, although the 
excitation mechanism is still uncertain and may involve IR 

In particular, see Fig. 4 of iMenten & Melnickl Il99lh . which 
shows that single tell-tale narrow (~ 1 km s~'-wide), long lived 
(> 10 y) feature at 22.3 km s _1 is seen in SiO maser lines from all 
three Si isotopes, as well as in both V2 = l lines. In fact, strikingly, the 
4 40 - 5 33 H 2 and the 29 SiO and 30 SiO v=0, 7=1-0 lines show only 
emission in this feature. 
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661 - 7 52 


293664.4" 
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1.3(0.1) 


~ [11,32] 


22(3) 


0.055 


~ 21 


660 - 753 


297439.1" 


3920 


<0.16 b 






< 0.029 c 




523 ~ 6l6 

y = 8-7 


336227.6" 


2939 


1.2(0.1) 


~ [11,32] 


18(2) 
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~ 21 


29 SiO v= 3 


335880.3 d 


5280 


0.78(0.07) 


[14,32] 


8.4(0.9) 


0.087 


23.5(0.3) 


30 SiO v= 1 


336602.4 d 


1804 


41.3(0.5) 


[-9,64] 


~ 2 C 


3.P 


19.5 C 


29 SiO v= 


342980.6 d 


58 


34.2(0.3) 


[-16, 62] 


39.6(0.5) 


0.81 


21.9(0.2) 



Columns are (from left to right) quantum numbers of upper and lower state, frequency, energy above ground of lower state, integrated main- 
beam brightness temperature, velocity range covered by line (FWZP), FWHM line width from a Gaussian fit, peak main-beam brightness 
temperature, and centroid LSR velocity, Multiply TmbCK) by 30 to g et flux density in Jy. Numbers in parentheses are formal errors derived 
from Gaussian fitting and/or noise level analysis, fchen et alJ fcOOOh b 3<x upper limit assuming a 20 km s 1 wide velocity range. c 3cr upper 
limit ( in a spectrum smoothed to 1.5 km s" 1 resolution). d Calculated from high precision molecular constants measured by Mollaaghababa et al. 
Il99ll) . e Value for strongest velocity component. 



line overlaps. Could we, on the other hand, expect to observe 
thermal emission from a region of a size similar to that of 
the SiO maser region in lines unaffected by any inversion 
mechanism? 

4.2. The nature of the water lines detected with APEX 

To investigate whether the observed emission in one or both 
of the H2O lines detected by us might be of thermal nature, 
we make the following considerations. Assume that the line's 
level populations were in Local Thermodynamic Equilibrium 
(LTE), meaning T ex = T mt = T v n, = 7\j n , i.e., that excitation, 
rotational, and vibrational temperatures are all equal to the 
kinetic temperature in the emitting region and that the lines 
were thermalized, i.e. their level populations were described 
by a Boltzmann distribution at that temperature. If we assume a 
uniform medium, the line's main-beam brightness temperature, 
Tmb, is given just by / x 7\i n (l - e~ T ), where the beam-filling 
factor, /, is given by (#s/#b) 2 - #s and 9b are the FWHM size of 
the Gaussian source and beam. 

We now estimate whether plausible assumptions on the 
emission region could produce optically thick emission. 

Around VY CMa, vibrationally (v=l and 2) SiO masers 
arise from a region of radius ~ 0."15 (Menten & Reid, 
in pr eparation). A ssuming a distance of 1.5 kpc for VY 
CMa jLada& Reidll978t) . this corresponds to 3 10 15 cm or 20 
stellar radii. Inside that region, SiO pumping models require a 
molecular hyd rogen den sity, w(H?), betwee n 10 9 and 10 10 cm" 3 
(lLockett & EhtzuJl992l:lBuiarrabail994T) . 

In the following we assume a constant 
density of 2 10 9 cirT 3 within 10 15 cm. This seems to be a 
plausible ave r age v alue between n(H2) = 7 10 7 cm" 3 , the value 
Danc hi et al 1 dl994l) find for the inner radius of the dust shell 
and the densi ty of 10 12 cm " 3 in the stellar photosphere; see, 
e.g., iReid & Mentenl Jl997h . Danchi et al. derive T = 1563 K 
at the inne r dust shell radius, while iLe Sidaner & Le Bertrel 
i 19961) and iMassev et alJ (|2006) derive a temperature of 750 
- 850 K for the dust from IR photometry. In the following we 



adopt 1000 K for the temperature. Assuming a homogeneous 
medium, we obtain an column density of 2 10 24 cm -2 . 

If we assume that the observed V2=l lines are thermalized 
we now can calculate their optical depths. For this we 
have used A-values and level energies from the HITRAN 
database ( Rothma n et alJ E003I and a value for t he par tition 
function, 2(1000 K), of 1194 from IChen etail (H)00). We 
use the upp er range of the [H7O/H7I abundance ratios of 
4 10" 4 that Gonzalez- Alfonso & Cernicharol ( 1 19991) derive for 
stars w ith high mass-loss rates while noting that Izubko et alJ 
(2004J) derive 4 10" 4 for VYCMa specifically in a completely 
independent analysis. 

For T = 1000 K we compute optical depths of 1 .0, 0.36, and 
2.8 for the 294, 297, and 336 GHz lines, respectively, and Tmb~ 
values of 12, 5, and 28 mK. The lower value of the 297 GHz 
para-line as compared to the 294 GHz ortho-line is due to the 
3:1 ortho-to-para ratio. 

The calculated value for the 336 GHz lines is only a 
factor of 2.3 lower than the actually measured intensity, while 
that ratio for the 293 GHz line is larger, ~ 5. For the 
broa d component of the 232 .6 GHz 550-643 line detected 
by iMenten & Melnick Jl989h with a 12" FWHM beam we 
calculate Tmb =43 mK for T = 1000K, which is smaller than 
the 140 mK they observe. As discussed in Sect. 14. H one might 
expect inversion for both of the latter lines and the weak maser 
action they do indeed seem to show (see Sect. |3j can explain 
their intensities. 

The above estimates assume that our H2O lines are 
thermalized. Thermalization can either be achieved by 
collisional excitation or radiative excitation. To check the 
potential of collisional excitation, we calculate the lines' 
critical densities. The critical density of a transition from level 
i to level j is given by £ Ay/ £ Jip where the sums are over all 
possible transitions originating from level i. Ay are the Einstein 
A-coefficients (again from HITRAN) and y,-,- t he collisional 
rate coefficients, calculated bv lGreen et al.l Jl993b for different 
T-values and access ed online via the Leiden Atomic and 
Molecular Database; ( Schoi er et al J2005I) . Rate coefficients for 
collisions from the ground to the V2 = 1 state are unknown. In the 
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following we assume them to be identical to transitions within 
the ground state with identical J, K a a nd K c quantum numbers 
multiplied by a factor c vg ; the study of iGonzalez-Alfonso et alJ 
(2002) suggests c vg =0.02. For all of our three V2=l transition 
we calculate 2A,y = 21s~ 1 within 10%. The sums of the 
collisional rate coefficients for the three lines are very similar as 
well: we calculate (within 5%) Eyo' = c vg 2.1 10~ 10 cm 3 s _1 for 
1000K. From this we derive re C rit = cVg 10 11 cirT 3 . Comparing 
this with the densities discussed above (~ 10 8 — 10 10 crrr 3 ) 
we find that our lines are very unlikely to be thermalized by 
collisions in the whole 10 15 cm radius region. 

As to radiative excitation: 

According to ICarroll & Goldsmith! J 1 98 lit , in a region with a 
strong IR field the criterion for radiative excitation to determine 
the populations of energy levels of a rotational line that is 
coupled by IR radiation of frequency vir to levels in an excited 
vibrational state is f/( e hym/kTs - 1) > A mt /A vib . Here A vib and 
A rot are the A-values the rotational and vibrational transitions 
involved, respectively, and T$ is the effective temperature of the 
IR source. For our case the filling factor, /, is equal to 1 and 
vir ~ 4.5 10 4 GHz (corresponding to 6.7 yum), which results in 
the ratio in left hand side of the equation to being « 0.2. This 
is much larger than the right hand side, since A V1D ~ 20 s _1 and 
A ro t ~ 10~ 6 s~'. Given this, it seems safe to assume that the 
level populations of the submillimeter lines are thermalized at 
the temperature of the emitting region. 

In summary, the level populations of the observed H2O 
V2=l lines are likely to be excited via infrared radiation. 
Making simple, but plausible assumptions on their dense, 
hot emission region, we can roughly reproduce the observed 
intensity of the 336 GHz line assuming thermal excitation. 
Known characteristics of H2O V2=l masers lead one to 
expect that the 294 and 297 GHz lines might have inverted 
populations. Indeed, the observed intensity of the 294 GHz 
ortho-line suggests that it might be boosted by weak maser 
action. The non-detection of its twin para-H20 297 GHz line is 
naturally explained by the lower abundance of the para species 
and less amplification. Future imaging with ALMA, whose 
highest resolution will resolve the lines' emission region, 
will yield important information on this higly interesting 
environment. 
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Fig. 1. Top to bottom: Spectra observed toward VY CMA of 
the H2O V2=l, 523 - 616, 661 - 752, and 660 - 753 lines and the 
7=8-7 lines of 29 SiO in the v=3 and 1 states and 30 SiO in the 
v=l state. The spectra have been smoothed to channel spacings 
of 1.3, 1.5, 1.5, 1.3, 1.3, and 0.22 km s , respectively. 
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